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Abstract 



We study some spring mass models for a structure having a unilateral spring of small 
rigidity e. We obtain and justify an asymptotic expansion with the method of strained 
' coordinates with new tools to handle such defects, including a non negligible cumulative 

(-H I effect over a long time: ~ e^^ as usual; or, for a new critical case, we can only expect: 

~ We check numerically these results and present a purely numerical algorithm 

to compute "Non linear Normal Modes" (NNM) ; this algorithm provides results close to the 
asymptotic expansions but enables to compute NNM even when e becomes larger. 
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^ ■ 1 Introduction 

O 

For spring mass models, the presence of a small piecewise linear rigidity can model a small 
defect which implies unilateral reactions on the structure. So, the nonlinear and piecewise 
. linear function u+ = max(0,u) plays a key role in this paper. For nondestructive testing we 

^ ' study a singular nonlinear effect for large time by asymptotic expansion of the vibrations. 

New features and comparisons with classical cases of smooth perturbations are given, for 
instance, with the classical Duffing equation: il + u + ev? = and the non classical case: 
ii + u + eu-|_ = 0. Indeed, piecewise linearity is singular: nonlinear and Lipschitz but not 
differentiable. We give some new results to validate such asymptotic expansions. Further- 
more, these tools are also valid for a more general non linearity. A nonlinear crack approach 
for elastic waves can be found in [llj . Another approach in the framework of non-smooth 
analysis can be found in [2l |4l [19] . 

For short time, a linearization procedure is enough to compute a good approximation. But 
for large time, nonlinear cumulative effects drastically alter the nature of the solution. We 
will consider the classical method of strained coordinates to compute asymptotic expansions. 
The idea goes further back to Stokes, who in 1847 calculated periodic solutions for a weakly 
nonlinear wave propagation problem, see [15l [161 HZl [H] for more details and references 
therein. Subsequent authors have generally referred to this as the Poincare method or the 
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Lindstedt method. It is a simple and efficient method which gives us approximate nonhnear 
normal modes with I or more degrees of freedom. 

In section [2] we present the method on an explicit case with an internal Lipschitz force. We 
focus on an equation with one degree of freedom with expansions valid for time of order e^^ 
or, more surprisingly, e^^/^ for a degenerate contact. 

Section [3] contains a tool to expand (u + and some accurate estimate for the remain- 
der. This is a new key point to validate the method of strained coordinates with unilateral 
contact. 

In Section m we extend previous results to systems with N degrees of freedom, first, with 
the same accuracy for approximate nonlinear normal modes, then, with less accuracy with 
all modes. We check numerically these results and present a purely numerical algorithm to 
compute "Non linear normal Modes" (NNM) in the sense of Rosenberg 2^ ; see PP for two 
methods for the computation of NNM; see [9] for a computation of non linear normal mode 
with unilateral contact and [14j for a synthesis on non linear normal modes; this algorithm 
provides results close to the asymptotic expansions but enables to compute NNM even when 
e becomes larger. 

In Section [51 we briefly explain why we only perform expansions with even periodic functions 
to compute the nonlinear frequency shift. 

Section [6] is an appendix containing some technical proofs and results. 

2 One degree of freedom 
2.1 Explicit angular frequency 

We consider a one degree of freedom spring- mass system (see figure [T]): one spring is classical 
linear and attached to the mass and to a rigid wall, the second is still linear attached to a 
rigid wall but has a unilateral contact with the mass; this is to be considered as a damaged 
spring. The force acting on the mass is kiu + /c2'"+ where u is the displacement of the 
mass, fci, the rigidity of the undamaged spring and k2, the rigidity of the damaged unilateral 
spring. We notice that the term 

= max(0,M). 

is Lipschitz but not differentiable. Assuming that fc2 — eki and dividing by the value of the 
mass we can consider the equation: 

u + cjgit + eu+ = 0, (1) 

where luq a positive constant. This case has a constant energy E, where 




Figure 1: Two springs, on the right it has only a unilateral contact. 

Therefore, the level sets of E{u, ii) will be made of two half ellipses. Indeed, for u < the 
level set is an half ellipse, and for w > is another half ellipse. Any solution u{t) is confined 
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to a closed level curve of E{u, ii) and is necessarily a periodic function of t. 
More precisely, a non trivial solution {E > 0) is on the half ellipse: ii^ + uj^u = 2E, in the 
phase plane during the time Tc = tt/ojoi and on the half ellipse + (luq + e)u = 2E during 
the time Te = n/ ^/uj^+e. The period P{e) is then 

— 1/2 

P{e) = (1 + (l + e/lUq) )n/ujo, and the exact angular frequency is: 
uis) - 2LUoil + {l + e/LU^,y'^Y^ 

Let us compare with the angular frequency for Duffing equation where the nonlinear term is 

u + WqU + eu^ = 0, (3) 
which depends on the amplitude uq of the solution ( see for example [15 1 ITS l [T7 l I18j): 

2.2 The method of strained coordinates 

Now, we compute, with the method of strained coordinates, LOg, an approximation of the 
exact angular frequency uj{e) which is smooth with respect to e by exact formula 

uj{e) = uj, + Oie^). 

We expose this case completely to use the same method of strained coordinates later when 
we will not have such an explicit formula. 

Let us define the new time and rewrite equation ([T]) with the new time 

S = UJet, u^{t) = Ve{s), (4) 
u;lv':{s) + u;lv,{s) + eiv,{s))+ = 0, (5) 

To simplify the exposition and the computations, we take following initial conditions for 

u,{0) = ao>0, iieiO) = 0, (6) 

i.e. We(0) — aQ and v'^{0) = 0. Similar computations are valid for negative oq, see Proposition 
l2.1l below. With more general data, i.e. when Ue(0) ^ 0, computations are more complicate 
and give the same angular frequency uj^, see section [5l 
In new time s, we use the following ansatz 

Veis) = uo(s) +ewi(s) + eVe(s). (7) 

and the following notations: 

uje = ujQ + euji + e^Lj2, uj'^ = ao + eai + e^a2 + 0{e^) (8) 
ao = Lul, ai ^ ^ivoUJi, a2 ^ ujf + 2ujouj2, (9) 

where tt'i,a;2 or ai,a2 are unknown. 

We will also use the natural expansion, {u + ev)+ = u+ + eH{u)v + • • • , where H is the 
Heaviside function, equal to 1 if u > and else 0. This expansion is validated in Lemma ISTTl 
below. 
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Now, replacing this ansatz in ([5]) we obtain differential equations and initial data for vo,vi,rg; 
set 



L{v) ^ -aoiv" + v) (10) 

Livo)= 0, vo{0) = ao,v'o{0) = 0, (11) 

L{vi)= (t;o)++ai<, fi(0) = 0, v[{0)^0, (12) 

L{r,)= H{vo)vi+a2v'^ + aiv'{ + R,{s), r,(0) = 0, <(0) = 0. (13) 

We now compute, ai, vi and then a2. We have vo(s) — aq cos(s). A key point in the method 
of strained coordinates is to keep bounded vi and for large time by a choice of ai for ui 
and a2 for r^. For this purpose, we avoid resonant or secular term in the right-hand-side of 
equations ([72]) . (I13p . Let us first focus on ai. Notice that, vo{s) = aocos(s) and aq > 0, so 



{vo)+ = ao 

I cos(s)| has no term with frequencies ±1, since there are only even frequencies. Thus ((uo)+ — 
ai?;o) = a-o cos(s)(l/2 — ai) -I- ao\ cos(s)|/2 has no secular term if and only if ai = 1/2, so 
uJi = l/(4(jJo)- Now, vi satisfies: 

-Luliv'l + vi) ^ ao\ cos s\/2, ui(0) = 0, v[{0) ^ 0. 

To remove secular term in the equation (llSp we have to obtain the Fourier expansion for 
H{vq) and vi. Some computations yield: 



2 4 ;^ (-i) fc 



^^0^(^)1 = ---Eifc^'=o«(2M, (14) 

fe=i 



.i(.) = _|[^l-25:^Xl)_eos(2Mj-^cos(.), (15) 

H{vo) = --f-^i^cos((2j--M).), (16) 

fe=i 



where A = ^ (1-2^ 



To remove secular term of order one in (fT^. it suffices to take a2 such that: 

[H{vQis))vi{s) + aawo (s) + ai<(s)] ■ wo(s)ds (17) 







For Duffing equation, see [151 HSl Hi] i the source term involves only few complex exponentials 
and the calculus of a2 is explicit. For general smooth source term, Fourier coefficients decay 
very fast. Here, we have an infinite set of frequencies for vi and H{vo), with only a small 
algebraic rate of decay for Fourier coefficients. So, numerical computations need to compute 
a large number of Fourier coefficients. For our first simple example, we can compute explicitly 
a2- After lengthy and tedious computations involving numerical series, we have from (jl7p 
and p^ . p5)) . to evaluate exactly the following numerical series: 



a2 = — 
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Figure 2: vi, logig ( — ) 
\max \Us\ J 



thus W2 = —{2ujo) as we have yet obtained in ([2]). 

When oq = 1 we obtain in figure [2] first modes of the infinite Fourier spectra for 

VoiuJet) + evi{ujet) ~ 

Indeed we have for negative or positive uq the following result. 
Proposition 2.1 Let be the solution of with the initial data: 

Ue{0) = oq + sai, Uj(0) = 0, 

then there exists 7 > 0, such that, for all t < — js~^ , we have the following expansion: 

Ue (t) = Wo {uJet) + evi {u;,t) +0{e'^), 



viiujet) = ( 1 - 2 T^p^;^— !— 2 cos(2fcwe0 ) + (ai - A)cos(wei), 



^hereA^^(l-2Y, 




In the Proposition 12.11 with the method of strained coordinates, we recover an asymptotic 
expansion for the exact angular frequency uj{e) = ljq + ujiE + ijJ2s'^ +0{e^) and for the exact 
solution Ue{t). The term vi is explicitly given by its Fourier expansion. Notice also that L02 
is not so easy to compute. It needs to compute a numerical series. 
The technical proof of the Proposition 12.11 is postponed to the appendix. 
Examples from Proposition l2.1l have angular frequency independent of the amplitude. Equa- 
tion ll]) is homogeneous . Indeed, it is a special case, as we can see in the non homogeneous 
following cases. In these cases, we assume that the spring is either not in contact with the 
mass at rest (6 > 0) or with a prestress at rest (6 < 0). 
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Proposition 2.2 (Nonlinear dependence of angular frequency ) 

Let b be a real number and let be the solution of: 

il + uJqU + ea{u — b)+ = 0, ^^(O) = ao + eai, Us{0) = 0. 

If \ao\ > \b\ then there exists 7 > 0, such that, for t G [OjT^], with Tg = —, we have the 
following expansion in C^([0,Te],R); 

Us{t) = vo{ij0^t) + evi{LList) + 0{e'^) with 

vq{s) = aocos(s), 

+00 

^^{^) = y^cifc cos(A:s), 
fe=i 



do = [sin(/3) — - — -p\ where (i — arccos [ — ) e [0,7r], 



dk 



-a\ao\ / siii((/c + l)/3) sin((fc - 26sin(A:/3) \ 
■KUjl{l-k-^)\ k + l k-1 \ao\k J' - ' 



di = ai — A where A = dk, 

We = LOo + eUJi + e'^UJ2, 

a /sin(2/3) ^ 2&sin(/3) 



2TTijjQ \ 2 |ao| 

W2 = ^^^^ — / H{aocos{s) — b)vi{s)cos{s)ds. 

ao woTrao Jq 

Notice that if |oo| < there is no interaction with the weak unilateral spring. Thus the 

linearized solution is the exact solution. 

Proof : There are two similar cases, ao positive or negative. 

First case : assume ao > 0. With the previous notations, the method of strained coordinates 
yields the following equations: 

Vq + uo = 0, wo(0) = ao, ?)o(0) =0 so va{s) = ao cos(s), 
—at){v'l + vi) = a{vQ — — cxiVo = aao(cos(s) — 6/ao)+ — aiao cos(s), 
-ao{r^ + rs) = aH{va - b)vi - a2Va ~ aivi + R^. 

After some computations the Fourier expansion of (cos(s) — c)+ for |c| < 1 is: 

+00 

(cos(s) — c)+ = Cfe cos(A;s), (18) 
fe=o 

sin(/3) — c/3 , ^ / ^ r« T 

Co = where p = arccos (c) e [0,7rJ, (19) 

TT 

CI = i(^+/3-2csin(/3)), (20) 
_ 1 / sin((fc + l)/3) , sin((fc-l)/3) 2csm{k(i) \ 

/■TT 

The non secular condition / {a{vo — 6)+ — aivo) cos{s)ds = 0, with c = b/ao yields ai = aci . 

Jo 

Now, we can compute coi = ai/(2a;o) and vi with a cosines expansion: vi{s) = ^ dk cos{ks) 
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with dk = for k ^ 1. The coefficient di is then obtained with the initial data 

ao 1 — fc^ 
i;i(0) = ai,{;i(0) = 0. 

a2, is obtained with the non secular condition for Tg: 

I r 

— — {aH{vQ — b)vi — a2VQ — aiWi) cos{s)ds. This condition is rewritten as follow 
I" Jo 

/ ij(aocos(s) — fe)wi(s) cos(s)(is, which gives uj2 since uj2 



a2 — 

ao 

Second case : when oq = — |ao| < 0, by a similar way, we obtain a similar expansion, except 
that (wo(s) — — |ao|(— cos(s) — 6/ao)+. And we only need the Fourier expansion of 
(— cos(s) — c)+ — \^ Cfe cos(fcs). 



k 



Co 



sin(/3) + c/3 



TT 



1 /sin(2/3) ^ „ . 
ci = — ^+/3 + 2csin(/3) 

TT \ Z 

1 ^ sin((fc + l)/3) , sin((fc~l)/3) , 2csin(fc/3) -, 



□ 



When |ao| = \b\, we have another asymptotic expansion only valid for time of order 
when the unilateral spring interacts with the mass. 

Proposition 2.3 (The critical case ) 

Let b be a real number, b ^ 0. and consider, the solution of: 

il + ujqU + ea{u — 6)+ — 0, Me(0) = ao + eai, 'ite(O) = 0. 
If 1^0 1 — \b\ then we have 

Ue(t) = (ao +£ai)cos(cjoi) +C(£^), 

— if |ao + eail > |6| where 7 > 0, 
+CX) eZse. 

The method of strained coordinates gives us the linear approximation for u^(t), with s ~ t, 
i.e. uje ~ 1- If l*^e(0)| < \b\, the exact solution is the solution of the linear problem H+ujqU — 0. 
Otherwise, if |we(0)| > since, |6| is the maximum of vo{s) — ao cos(s), a new phenomenon 
appears, during each period, |ue(t)| > |6| on interval of time of order ^/s instead of e. Then 

is smaller than in Proposition [2711 
To explain this phenomenon, we give precise estimates of the remainder when we expand 
{vq + evi + e'^re)+ in the next section, see Lemmas l3. II and l3^ below. 



3 Expansion of {u + ev] 



We give some useful lemmas to perform asymptotic expansions and to estimate precisely the 
remainder for the piecewise linear map u — > m+ = max(0, u). 

Lemma 3.1 [Asymptotic expansion for (u + £u)+ ] 
Let be T > 0, M > 0, u,v two real valued functions defined on I = [0, T], 

Je = {te I, \u{t)\ < eM}, 

tJ-e{T) the measure of the set and H be the Heaviside step function, then 



{u + ev)+ = {u)+ + sH(u)v + exe{u,v), with H(u) 



1 ifu>Q 
elsewhere 
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and Xe{u,v) is a non negative piecewise linear function and 1-Lipschitz with respect to v, 
which satisfies for all e, 
'if\v{t)\ < M for any t G I: 



T 



Xe{u,v)\<\v\<M, / \Xe{u{t),v{t))\dt<M^Xe{T). (22) 



The point in inequality ()22|1 is the remainder exe is only of order e in L°° but of order e/ig 
in . In general, ji^ is not better than a constant, take for instance m = 0. Fortunately, it 
is proved below that /i^ is often of order e, and for some critical cases of order ^/e. 
Proof : Equality (|22p defines Xe and can be rewritten as follow: 

X.{u,v) = {u + ev)^-u^-eH{u)v _ ^^3) 

So, Xe is non negative since u — > u+ is a convex function. We also easily see that the map 
{u, v) — > Xs{u, v) is piecewise linear, continuous except on the line u — where Xe has a jump 
—V. This jump comes from the Heaviside step function. An explicit computations gives us 
the simple and useful formula: 

/ \ ( \u + ev\ if \u + ev\ < \ev\ 

= I elsewhere " ^^^^ 

We then have immediately < Xe{u,v) < \v\. Let u be fixed, then v Xe{u,v) is one 
Lipschitz with respect to v. Furthermore, the support of Xe is included in J^, which concludes 
the proof. □ 

We now investigate the size of /i^ (T) , see [S] [10] for similar results about /i^ (T) and other 
applications. 

With notations from Lemma l3JJ we have. 



Lemma 3.2 (Order of pL^{T)) Let u be a smooth periodic function, M he a positive con- 
stant and ^e{T) the measure of the set {t £ /, \ u{t)\ < eM}. 
If u has only simple roots on I — [0,T] then for some positive C, 

fie{T) < CexT. 

More generally, if u has also double roots then 

tie{T) < CV^xT. 

The measure of such set Jg implies many applications in averaging lemmas, for a character- 
ization of fie in a multidimensional framework see [3l [10] . 

Notice that any non zero solution u(t) of any linear homogeneous second order ordinary 
differential equation has always simple zeros, thus for any constant c the map t ^ u(t) — c 
has at most double roots. 

Proof : First assume u only has simple roots on a period [0, P], and let Z = {to G 
[0, P],u(to) — 0}. The set Z is discrete since u has only simple roots which implies that 
roots of u are isolated. Thus Z is a finite subset of [0, P]: Z = {ii, t2, • • • , ijv}- We can 
choose an open neighborhood Vj of each tj such that w is a diffeomorphism on Vj with 
derivative |u| > \u{tj)\/2. On the compact set K = [0, P] — UV,, u never vanishes, then 

AeM 

min \u{t)\ — £q > Q. Thus, we have for all eM < Eq, the length of in is | V^n J^j < — — —. 
teK \u[tj)\ 

As jie is additive {fieiP + t)^ /ie(P) + l^eit)), its growth is linear. Thus, for the case with 

simple roots, we get fJ-dT) — 0{eT). 

For the general case with double roots, on each small neighborhood of tj-. Vj, we have with 
a Taylor expansion, |u(ij + s)\ > dj\s\'- , with I <l <2, dj > so, |V, n Je| < 2{eM/djY/\ 
then lJ,e{P) = ©(a/e), which is enough to conclude the proof. □ 
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4 Several degrees of freedom 



Now, we investigate the case with N mass. We use, the method of strained coordinates in 
three cases. We present the formal computations for each expansion. The mathematical 
proofs are postponed in the Appendix. 

In subsection 14.11 the initial condition is near an eigenvector such that the approximate 
solution stays periodic. We gives such initial condition near an eigenvector in subsection 14.21 
to get an approximate nonlinear normal mode up to the order e^. Finally, in subsection 14.41 
all modes are excited. An extension of the method of strained coordinates is still possible 
but only at the first order with less accuracy. 
The system studied is the following: 

MU + KU + e{AU — = 0, where, for each component, 

/N \ 

[{AU — B)^]k = ^^OfcjUj — &/c , M is a diagonal N x N mass matrix with positive 

\^=i 7 + 

terms on the diagonal, K is the stiffness matrix which is symmetric definite positive. It is 
also possible to add many terms e{AU — B)+ modeling small defects. For a such system, 
endowed with a natural convex energy for the linearized part, we can control the e-Lipschitz 
last term for e small enough up to large time, so for e << 1 the solutions remain bounded 
for time of the order . 

To simplify the exposition, we deal with the following diagonalized system for the linear part, 
keeping the same notations, except for the positive diagonal matrix A: 

il + K^U = -e{AU-B)+. (25) 

In fact, with a linear change of variables, we have in the right hand side of the previous 
equation a linear combinations of terms £{AU — -B)+. 

4.1 Initial condition near an eigenvector, 
second order approximation 

For the system (j25p . we take an initial condition near an eigenmode of the linearized system 
denoted for instance by index 1 . 

uf(0) = ao + eai, uf (0) = 0, 
u|(0) = +eafe, u|(0) = 0, fork^l. 

We impose 02 , • • • , a^r later to have a periodic approximation, but ai is a free constant as 

cq. It is a key point to apply the method of strained coordinates. 

We use the same time s — ojet for each component and the following notations. 

L^J^ = Wo + £Wl + £^^2, Wo = Ai, 

{uJeY = ao + eai +£2a2 + 0(£^), ao = wg = Af, 

ai — 2woWi, a2 = (jjf + 2ujoUJ2, 

up) = vl{s)^v^is) + ev]is) + e^r'jis), j = l,---,iV. 

Replacing, this ansatz in the System (|25p we have in variable s, 

N 
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and then performing the expansion for all fc G {1, • • • , N}, 
First we have Vi{s) = aq cos(s). 

Equations for v'^, for all k ^ 1, with zero initial data give us w° = 0. 
In equation for vl, we remove the secular term in the right hand side, 

-aoiivlr + vl) = (anf? - bi) ^ + = r.h.s. vl{Q) = a,, K)'(O) = 0. 

The orthogonality of the r.h.s with cos(s) leads to define ai. For instance, if bi — 0, we have 

as in Proposition [221 Q^i = -rr — 7T~- 

Now, ai is fixed, so vl is a well defined even 27r periodic function. 
Then, for fc 7^ 1, we can compute with the simplified equation, 

ao{vl)" + Xlvl = -(afci«?-6,.)+- (26) 

Let be the unique 27r periodic solution of equation (|26p . Such function exists and is 
unique if Xk ^ AiZ. Furthermore is an even function as the right hand side of equation 
(EH). Let vl be (111, i.e. vl{0) = ak = 0^(0) and (t;^)'(0) = for all k ^ I. 
The term rf, with null initial data, has a simplified equation since v^ = for all A; 7^ 1, 

-Lirf = if (ani'J-foi) |^E«ij"lj +"2K)" + aiK)" + i?!. 

We now can compute numerically 02 to avoid secular term in the right hand side, i?f ex- 
cepted, with the following condition. 



N 



■ cos(s)ds. 



\k=l / 

Rewriting this condition, we obtain an equation for a2 in Thcorem l4.1l below. which gives 
UJ2 as in the Proposition [^21 

For each fc 7^ 1, Afe and Ai are Z independent. So r| stays bounded for large time since there 
is no resonance of the order one at the first order in equations, —Lkrf. = (•••) + -^1- This 
is the technical part of the proof to validate rigorously and to find the time of validity of 
such asymptotic expansion. The complete proof to bound (rf, • • • ,r^) for large time is to 
be found in the Appendix. 

We now state our result with previous notations. 

Theorem 4.1 The following expansion is valid on {O^T,,), with — > +00 when e ^ 
under assumption Afc and Ai are Z independent for each k ^ 1, 

u\{t) = vl{LO,t) + ev\{uo,t) + 0{e^), , . 

u%it) = + eviico.t) + k^l, ^^'> 
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where Vi{s), ai, uji, v\{s), ^^(s), 4>\, a2, ^2 o-fe defined successively as follows, with differ- 
ential operators = ^I'J^ + -^L 

v^{s) = aocos(s), (28) 

ai = / {aiiVi{s) — bi)+ cos{s)ds, (29) 

- - (^«) 

-L^vl = (ant'?-&i)+ + ai(«?)", vHO) ^ a,, {vD' (0) = 0, (31) 

-Lkvl = (afei«?-6fe)^, «i(0) = a, :=0i(O), («i)'(0) = 0, /or/c^l, (32) 

is t/ie unique 2n periodic solution of \3l3jl for all k ^ 1, (33) 



2 r \ I ^ 

— { H{aiiv1{s) - bi) aiivl + aut/)^ 



= I ~y <( iJ(aii?;i(s) - 61) I aiiui + 2^aifc(/)fc I cos(s)ds I - ai, (34) 

UJ2 = • (35) 

Furthermore, if (ajiWi — bj) has got only simple roots for all j = 1, ■ ■ ■ , N , 
then is of the order 
else is of the order 

In the theorem is classically obtained by a Fourier series. We give some indication of its 
initial condition in the next subsection 14.21 

4.2 Approximate non linear normal mode 

The special initial conditions of the previous subsection can be explicited in order to find 
a solution where all the components are in phase at the same frequency. Indeed we shall 
obtain an approximate curve of initial conditions for which the solution is periodic up to the 

order e for a time of the order — : this is up to the approximation a non linear normal mode 

e 

in the sense of Rosenberg see [S] for a computation of non linear normal mode with 
unilateral contact and [14] for a synthesis on non linear normal modes. 

Theorem 4.2 For the system I125\) . we take an initial condition close to an eigenmode of 
the linearized system denoted for instance by index 1 . 

■uf(O) = OQ + eai, Mf(0) = 0, 
u|(0) = O + eofc, u|(0)^0, fork^l. 

7 

the results of the previous theorem are still valid for ^ — but now we are looking for 

e 



ak = vl{Q) = in ([32|1. 

Moreover for the cases detailed below, the solution is periodic with angular frequency uj^ up 
to the order e and for < t < 

1. assume bk — 0, for k ^ 1 with the initial condition defined with 

_ Qfciao _ |QfciQo| afciQQ (-1)' /„„^ 

'''-2{Xl-Xl) Xln + 2 A.(4PAf-A2,)(4P-l) ^''^ 

2. assume < |^^''^'°^^| < 1, and aoOki < for k =/= 1 then for the initial condition defined 
with 



Ok — aoOki 



, +00 

Ofe Q 



(37) 
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where ci are defined in (11811 with c = — — 

J 1 ' a-kio-a 



3. assume < 
with 



\bk\ 
\a.kiao\ 



< 1, and aoflfei > for k ^ 1 then for the initial condition defined 



Qk = aoOki 



+ 00 

E 

.1=1 



Cl 



PA? 



(38) 



where ci are defined in psp with 



Remark 4.1 The other cases are less interesting but may solved similarly. 

Proof : The principle of the proof is simple; vl is solution of the differential equation 
with 'y^(O) = Ofc and has to be determined in order that the function vj. has an angular 
frequency equal to one. It is elementary that the solution of ([5^ is 



vj. = A cos 



B sin ( — 

,Al 



(39) 



where is a particular solution associated to the right hand side which is of angular fre- 
quency equal to 1; note that B = as the initial velocity is null; we can get a function of 
angular frequency equal to 1 by setting Ok = (f'kiO)', this condition may be written explicitly 
with formulas pS]) which provides the expansion in Fourier series; formulas ([36]), ([37|) . ([38]) 
are then derived easily. 

1. for 6fc = for A: 7^ 1, ([5^ is written: 



cos{s) 



\akiao\ 



. cos{s) 



we use formula (fT4|) to get the particular solution 



,1 cos{s) 
Ok = Ofciao 



2(A?-A^) A^ 



+ 00 

E 



(-1)' cos{2ls) 
2 ^^4PXl-Xll 4P-1 



(40) 



(41) 



from which ((36|) is deduced. 



2. for < 



-Lkvl 



|afciao| 



< 1, and aQOki < for k =/= I 



is written: 



-Lkvl = -afciflo [~cos{s) H ) 



-cos(s) 



-cos(s) 



cos(s) 



Ofciao 



-l-oo 



2(Af - A^) \%akiao 
h 



E 



cicos{ls) 

■2T2 



A 



(42) 

we use ([22]) to obtain 
(43) 

where q is defined with 

(44) 



(3 = arccos{- 



3. for < 



Idfciool 



< 1, and aoflfci > for fc 7^ 1 



) which defines c/ in ([22]) with c : 



Ofeiao 



-LfcWfc = flfciao ( cos(s) ) 



is written: 

from which 



+00 

E 



cicoslls) 



2 where c; is defined with 

I Ai — \i 



P = arccos(- 



which defines c/ in ([22]) with c : 



(45) 

(46) 
(47) 
(48) 
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□ 



4.3 Numerical results of NNM 

4.3.1 Using numerically Lindstedt-Poincare expansions 

Here we use the previous resuhs and compute numericaUy a solution of system (|25p using 
the approximation (|27p 

with the initial conditions of theorem l4.2l The first term v'^ is easy to obtain; for the second 
term an explicit formula is in principle possible using Fourier series such as for one degree 
of freedom but it is cumbersome so we choose to compute by solving numerically (j32p with 
a step by step algorithm; we use as a black-box the routine ODE of SCILAB ,25, to solve 
equations of theorem 14.11 after computing by numerical integration ai. We show numerical 
results for a system of the type: 

MX + KX + eF{X) = (49) 

we still denote the eigenvalues and the eigenvectors of the usual generalized eigenvalue 
problem 

K(t)j - \^jM(t)j = with (50) 
= 4, (51) 

We set: 

X = Uj = (f)u (52) 
j 

In this basis, the system may be written componentwise: 

Uk + Xluk +* 0fcF(£u) = (53) 

For a local non linearity in the system (|49p . written in the basis of the eigenvectors, we do 
not obtain a system precisely of the form (j25p : so we illustrate it with the non linearity: 

Fix) = MMXl - I3i) = Af0i(5] u,4>ji - (54) 

3 

where we denote by <j)i any normalized eigenvector of 1501 then the system (|49p is written: 

Uk + \\uk + - (55) 

k 

so that A and B of ([25]) are: 

ak] = bk = 

We find in figure [3] a numerical example of the Linstedt-Poincare approximation for 5 
degrees of freedom with e — 0.063 and with an energy of 0.03002. The left figure shows 
the 5 components of the solution with respect to time; the right figure, the solution in the 
configuration space: absissa component 1 and ordinate components 2 to 5; these lines are 
rectilinear like in the linear case but the non symmetry may be particularly noticed on the 
smallest component which corresponds to the mode where the non linearity is active. 
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Lind-Poincare, eps=0.063, Tmaxeps=3.2717 



0.06 
0.04 
0.02 
0.00 
-0.02 
-0.04 
-0.06 




Lind-Poincare eps=0.063 Tmaxeps=3.2717 vp1=3.6825 energie=0,03002 
0.06- 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 




-0.03 -0.02 -0.01 



,01 0.02 0.03 



Figure 3: Lindstedt-Poincare, energy=0.03, 5 dof; left: components with respect to time; ri 
in configuration space 



4.3.2 Using optimization routines 



We find also in figure 14.3.21 a numerical example with the same energy of 0.03002; it is 
computed with a purely numerical method described below. We notice that the solution is 
quite similar in both cases. 

The numerical expansions of the previous subsection gives valid results for e small enough; 
in many practical cases such as [7], e may be quite large; in this case, it is natural to try 
to solve numerically the following equations with respect to the period T and the initial 
condition X{Q). 



X{0)=X{T), X{0)^X(T) E{X)^e 



(56) 



In other words, we look for a periodic solution of prescribed energy; this last condition is 
to ensure to obtain an isolated local solution: the previous expansions show that in general, 
the period of the solution depends on its amplitude prescribed here by its energy. To try 
to solve these equations with a black-box routine for nonlinear equations such as "fsolve" 
routine of SCILAB [5S] (an implementation of a modification of Powell hybrid method which 
goes back to [5D]) in general fails to converge. Even in case of convergence, we should address 
the question of link of this solution with normal modes of the linearized system. 

So we prescribe that e — ce and for e — > 0, the solution is tangent to a linear eigenmode. 
In the case where all the eigenvalues of the linear system are simple, we define N (the number 
of degrees of freedom) non linear normal modes for which, it is reasonable to conjecture that 
they correspond to isolated solutions of (|56l) at least for small e if we enforce for example 
X(0) = 0. 

Algorithm This definition of the solution of tangent to a prescribed linear eigenmode 
provides a simple way of numerical approximation: using a continuation method coupled with 
a routine for solving a system of non linear equations. Define: 



T{e, Xo, T) = [X{T) - Xo, E{X) - ce] 
where X is a numerical solution of the differential system 
MX + KX + eF{X) = 
X{Q) = Xo, X{0) = Xi 



(57) 
(58) 
(59) 
(60) 
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Figure 4: Continuation and Powell hybrid, energy=0.03, 5 dof; left:fourier transform ; right: in 
configuration space 



choose a small initial value of e and an increment S 

choose an eigenvector 

Xo(0) = A,(t>,, Xi(0) = B,Xj(t>, 

with E{Xq{0),Xi{0)) = ce 

for iter=l:itcrmax 

e = £ + S 

with Xo(iter — 1) as a first approximation, solve for Xo(zter), 

^(e,Xo,r) = (61) 

if \\T{e,Xo,T)\\>tolerancethene^e-5, 5 = 5/2 
endif 
endfor 

This algorithm may be improved by using not only the solution associated to the previous 
value of £ to solve 

T{£,Xo,T) = Q (62) 
but also the derivative of the solution with respect to Xg, T. 

Numerical results These results are obtained by solving the differential equation with 
a step by step numerical approximation of the routine ode of Scilab without prescribing the 
algorithm. As we are looking for a periodic solution, this numerical approximation may be 
certainly improved in precision and computing time by using an harmonic balance algorithm. 
In figureOEl the same example with 5 degrees of freedom and energy equal 0.123 and 0.201 
are displayed. 

On the left of figure [5] we find the decimal logarithm of the absolute value of the Fourier 
transform of the solution; the Fourier transform is computed with the fast fourier transform 
with the routine f ft of Scilab; we notice the frequency zero due to the non symmetry of the 
solution and multiples of the basic frequency; no other frequency appears; on the right the 
five components are plotted with respect to time; we still notice the non symmetry. 
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On the left of figure [S] we find the solution in the configuration space and on the riglit 
the five components are plotted with respect to time; we still notice the non symmetry. 

In figure[7]we find results with 20 degrees of freedom, e — 0.272 and energy of 0.129; the 
NNM is computed by starting with an eigenvector associated to the largest eigenvalue . We 
see on the left in the configuration space that the components are in phase and on the right, 
the Fourier transform shows zero frequencies and multiple of the basic frequency. 

esp=0.06 epsu0=0.2593 Tper=3.2685 Tmax=65.3703 

15 T 1 

eps=0.06 epsu0=0.2593 Iam1=3.6825 ener=0.1234 Tper=3.2685 Tmax=65.37 




'O.O 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



Figure 5: energy=0.123, 5 dof; left:fourier transform ; right: vi^ith respect to time 




Figure 6: energy=0.2, 5 dof; left: configuration space; right:with respect to time 

In figure [8] the energy is 0.29 and the NNM is computed by starting with an eigenvector 
associated to the smallest eigenvalue; we notice on the left, the solution in the configuration 
space: at zero each dof has a discontinuity in slope which is clear. 

In figure [H the shape of the NNM is displayed on the left for the NNM starting from 
the eigenvector associated to the smallest eigenvalue and on the right for the NNM starting 
from the second smallest eigenvalue. We notice that the shape is quite similar to the shape 
of the linear mode. 

In figure [ini the NNM is computed by starting at an eigenvector associated to the second 
smalest eigenvalue. We notice that the first NNM is more asymmetric than the second one; 
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Figure 7: energy=0.129, 20 dof; left:in configuration space; right: fourier transform 




Figure 8: energy=0.29, 5 dof; left: configuration space; right:fft 



4.4 First order asymptotic expansion 

In this subsection, we do not particularize the initial data on one eigenmode. We adapt the 
method of strained coordinates since all modes are excited. We loose one order of accuracy 
compared to previous results since each mode does not stay periodic and becomes almost- 
periodic. 

More precisely, the method of strained coordinates is used for each normal component, with 
the following initial data 

4(0)=afc, it|(0) = 0, k = l,---,N. 

Let us define N new times Sk — Ap and the following ansatz, 

A| = + eA^, A^ = Afe, 

Kit) = vliXlt) - vlisk) - v"^{sk)+erlisk). 
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mode shape 

esp=0.06 epsu0=0.688 Iam1=0.0058 energie=0.29136 



0.0 
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mode shape 
,esp=0.06 epsuO= 
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0.4 
0.2 
0.0 
-0.2 
-0.4 
-0.6 
-0.8 



Iam2=0.0526energie=0.28173 




2 4 6 



10 12 14 16 18 20 



Figure 9: 20 dof; left:energy 0.29 mode 1; right: energy 0.28, mode 2 




Figure 10: energy=0.28, 5 dof; left: configuration space; right:fft 



The function are easily obtained by the linearized equation. Indeed, the only measured 
nonlinear effect for large time is given by To obtain these TV unknowns, we replace 

the previous ansatz in the system (|25p . 



(A|)'(f|)"(sfc) + \lvk{sk) = -e ( E "^-^-^l ( Al^'^ j - ^'^ 

The right hand side is written in variable Sk instead of Sj. Performing the expansion with 
respect to epsilon powers yields 

Lkvl = {\lf{vl)"{sk) + Xlvl{sk)^Q, (63) 
-Ltrlisk) = (jl<^k,v'',[^^s^^h^ +2\k\l{vlr + RI. (64) 
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Noting that replacing Vj(sj) by Vj (yn-Sfc) in (|64p implies a secular term of the order et, 

since Sj = -r^Sk + 0{et), the functions Vj are smooth and the map S — > 5+ is one-Lipschitz. 
These new kind of errors 0{et) are contained in the remainder of each right hand side: 



N 



^ E(^l)'- (65) 

\ fc=i 



RUt) = Oiet) + Oie\r^\), 

If &fc = 0, we identify the secular term with the Lemma 16.51 below and the relation 5*+ = 
S/2 + 15*1/2. Then, we remove the resonant term in the source term for the remainder r|, 

which gives us Xl = -tt— • 

If bk 7^ 0, we compute Aj, numerically with the following orthogonality condition to cos(s) 
written in the framework of almost periodic functions, 




■ cos{s)ds. 



The accuracy of the asymptotic expansion depends on the behavior of the solution </> = 
(01,- •• jiPn) of the N following decoupled linear equations with right coefficients to 
avoid resonance 



LkMsk) = \y.<^kjv° {-^sk\ ~bk\ +2XkXl{vt)". (66) 




Furthermore each function r| depends on all times sj, j — 1, • • • ,N and becomes almost- 
periodic, i.e. r| — f|(si, • • • ,sn)- Thus the method of strained coordinates, only working 
for periodic functions, fails to be continued. 

Nevertheless, we obtain the following result proved in the Appendix. 
Theorem 4.3 (All modes) 

If Xl, - ■ ■ , Aat are Z independent, then, for any — o{e^^), i.e. such that 

lim = +00. and lim e x Tg = 0, 

we have for all k = 1, ■ ■ ■ , N , 

\im \\u%it)-v° (A|i)lk^,~(o,T.) =0 

where X'j. = Xk + ^A^, w^(s) — a^ cos(s), and Aj. is defined by: 

xl = TTT— / I y^afci^j ( W^fc 1 "^fe I cos(s)ds. 




2Afcao T— >-i-oo T 7o 



Notice that accuracy and large time are weaker than these obtained in Theorem 14.11 It is 
due to the inevitable accumulation of the spectrum near the resonance and the various times 
using in the expansion. On the other side we have the following direct improvement from 
the Theorem 14. II 



Remark 4.2 (Polarisation) // only one mode are excited, for instance the number 1, i.e. 
ai 7^ 0, ak = for all k ^ 1, then we have the estimate for all t G [0, e"^]: 

ul{t)= i;?(Afi) +0{e). 

ul{t) = +0(£) for all k ^ I. 
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5 Expansions with even periodic functions 



Fourier expansion involving only cosines are used throughout this paper. There is never 
sinus. In this short section we explain why it is simple to work with even periodic functions 
and we give some hints to work with more general initial data. 

First, we want to work only with co-sinus to avoid two secular terms. If we return to 
equation (|12p : —ao{vi + vi) = (wo)+ + ckiWg- A priory, we have two secular terms in the 
right hand side, one with cos(s) and another with sin(s). Only one parameter ai seems not 
enough to delete all secular terms. 

Otherwise, if vq gM., u, S are 27r periodic even functions, g G C°(M, M) such that 

= e^^{S{s) + g{u{s)))ds then the solution of 
Jo 

v'' + v = S{s)+g{u), -y(O) = -^0, -^'(0) = 0, 

is necessarily a 27r periodic even function. Since we only work with 27r periodic even functions 
we have always at most one secular term proportional to cos(s). 

We now investigate the case involving not necessarily even periodic functions. In general, 
Uq and is the solution of 

lie + Ue + sf{Ue) = 0, Ue(0) = Uq, Ue{0) = Uq- 

By the energy 2E — ii'^ + + eF{u), where F' — 2/ and F{0) = 0, we know that is 
periodic for e small enough, for instance with an implicit function theorem see |28j also valid 
for Lipschitz function [4] in our case. Denote by Tg the first time such that U£{t) = 0. Such 
time exists thanks to the periodicity of u^. Now, let Ue defined by Ue{t) = Ue{t + t^). Ue is 
the solution of 

Ue + U, + ef{U,) = 0, C/e(0) = = Ue(Te), Ue{0) = 0- 

The initial data Uq depends on the initial position and initial velocity of Ug through the 
energy, (f/g)^ + eF{UQ) = (uq)^ + (^o)^ + ^-^(^o)- ^'^^ instance, if Uq and m° are positive 
then Uq is positive and 

US = + + e(i^K) - F + (4)2^ + 0{e^) . 

We can apply the method of strained coordinates for U^ only with even periodic functions: 
Ue{t) — vo^LUet) + evi{uJet) + 0{e'^). The expansion obtained for by U^, with (j)^ — —uj^Te 
is: 

Me(t) = VQiujet + (f>e) + evi{uet + (f)e) + 0{e^), 

which is a good ansatz in general for u^, where vq and vi are even 27r— periodic functions. The 
method of strained coordinates becomes to find the following unknowns 0o, uji,(j3i, LU2,(j)2 
such that 

oje = ujQ + euji + e^cj2 + ■ • • , 

(j)e = (j)o + £01 + £^02 H . 

Indeed, we have two parameters to delete two secular terms at each step. If one is only 
interested by the nonlinear frequency shift, it is simpler to work only with cosines. 

Otherwise, if / is an odd function, we can work only with odd periodic function. It is 
often the case in literature when occurs a cubic non-linearity. See for instance [TSl HZl HH] 
for the Duffing equation, the Rayleigh equation or the Korteweg-de Vries equation. 
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6 Appendix: technical proofs 



We give some useful results about energy estimates and almost periodic functions in subsec- 
tion [6lT] Next we complete the proofs for each previous asymptotic expansions in subsection 
16.21 The point is to bound the remainder for large time in each expansion. 



6.1 Useful lemmas 

The following Lemma is useful to prove an expansion for large time with non smooth non- 
linearity. 

Lemma 6.3 [Bounds for large time ] 

Let We be a solution of 

w'e + We = S{s) + /e(.s) + ege{s,We), 

WeiO) - 0, w',{0) = 0. ^ ' 

If source terms satisfy the following conditions where A/ > 0, C > are fixed constants : 

1. S{s) is a 2TT-periodic function orthogonal to e*"*, and |>5'(s)| < M for all s, 

2. \fe\<M and for allT, I \fe{s)\ds < CeT {resp. Cy/eT), 

Jo 

3. for all R > 0: Mr = sup \ge{s,u)\ < oo, 

e6(0,l):S>0,-R>|M| 

that is to say that ge{s,u) is locally bounded with respect to u 
independently from e G (0, 1) and s G (0, +oo), 

then, there exists Sq > and 7 > such that, for < s < Eq, w^ is uniformly bounded in 

W^^°° (0, Te), where = - {resp. ^). 

e 

Notice that and g^ are not necessarily continuous. Indeed this a case for our asymptotic 
expansion, see Lemme TS.ll and its applications throughout the paper. But in previous sections 
the right hand side is globally continuous, i.e. S + f^ + sge{.,We) is continuous, so, in this 
case. We is C^. 



Proof of the Lemma 16.31 First we remove the non resonant periodic source term which is 
independent of e. Second, we get bound for Wg and w'^ with an energy estimate. Third, 
with equation (p7)) . we get an uniform estimate for in L°°{0, T^) and the W'^'°° regularity. 
Step 1: remove S 

It suffices to write w^ = wi + where wi solves the linear problem: 

w'l + wi=S{s), ^1(0) = 0,^(0) = 0. (68) 

wi and w[ are uniformly bounded in L°^'{0, +00) since there is no resonance. 

More precisely, wi = F{s) + Acos(s) + i?cos(s), where F is 2tt periodic. F is obtained by 

Fourier expansion without harmonic n ~ ±1 since S is never resonant: 

^-^ 1 — ^ 

F is uniformly bounded, with Cauchy-Schwartz inequality set Cg = — we 

obtain: 

II^'IIl- < ^ - C'o||'5'||L2(0,27r) < C'o||'5'||L~(0,27r)- 
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Similarly, set = r?\r? - IT^, we have \\F'\\l^ < Doll5'|lL~(o,27r)- 

Furthermore, = Wi{0) = F{0) + A, and = (wi)'(O) = F'{0) + B, then, A and B are weh 
defined. w{ is also bounded, i.e. there exists Mi > such that || wi || vyi ^Co.+oo) < Mi. 
Notice that from equation (|68|) . wi belongs in 

Then we get an equation similar to (|67p for with 5 = and the same assumption for the 
same and the new g^'. 5e(s, w) — ge{s, wi + w). 

{wD" + {wD = /,(s) +eg,(s,u;i), , . 

{wim^O, K)'(0)=0. ^^^^ 

Step 2: energy estimate 

Second, we get an energy estimate for w^. We fix i? > such that R is greater than the 
uniform bound Mi obtained for wl, Mi + p with p > 0. Let us define 

2Eis) = (K)'(,s))2 + iw^,)is)', m = sup E{r), 

0<T<S 

and Te be the first time T > such that 2E{T) > p^, i.e. p estimates the size of (wf) and 
("'I)'- 

Multiplying the differential equation by (1^2)'' have for all s < T < T^{p) the 

following inequalities since sup \{w2y {t)\ < y 2E{s), and / |/e(s)|(is < CeT, 

0<T<s * Jo 

E{s) = r f,{r){w^)\T)dT + e rg,{T,{wl){T)){w-)'{T)dT, 
Jo Jo 



< Cesd2E{s) + esMR'J2E{s), 



E{T) < CeT\l2E{T)+eTMR\l2E{T), 
eT > 



E{T)/2 



Mr + C 

Notice that if 2'E{T) < p^ for aU T > then = +00. The critical case is when is 
finite and E{T) approaches when T goes to T^{p). Thus we have > ^ 



2e{MR + C) 

and E{t) < 4 for t < T, = 1 with 7 = ^(IT^- 

The proof is similar when / \fe{T)\dT < C^feT then > — / -. □ 

Jo 2^j£[y/£MR + C) 



For completeness, we state a similar and straightforward version of Lemma 16.31 useful for 
systems. 

Lemma 6.4 [Bounds for large time for systems ] 

Let Ws = (uif , • • • , w^) he the solution of the following system: 

{Xifiwir + {Xkfwi = Su{s) + + eg%{s- w,), 
wm = o, K)'(0) = 0, fc = l,...,iV. ^'"^ 

If source terms satisfy the following conditions where AI > 0, C > are fixed constants : 
1. non resonance conditions with Sk{s) are 2'K-periodic functions and \Sk{s)\ < M, 

(a) 5*1(5) is orthogonal to e*'^, i.e. / 5i(s)e***ds = 0, 

Jo 

(b) Afc,Ai are Z independent for all k ^ 1, 
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2. \fl\ < M and for all T, [ \fe{s)\ds < CeT or Cy/eT, 

Jo 

3. for all R > 0; Mr = max sup \gl{s;u)\ < oo, 

then, there exists Sq > and 7 > such that, for < s < Eq, is uniformly bounded in 

(0,T,), where Te = - or 

e y^e 

Proof : First we remove source terms Sk independent of e setting w| = Wk,i + w| 2 where 
Wk,i is the sohition of 

Xjwl, + Xlwk,i = Sk, Wfe4(0)=0, <i(0) = 0. 

As in the proof of Lemma [^751 wi^i belongs in W^'°° thanks to the non-resonance condition 
l.(a). For fc ^ 1, there is no resonance since ^ Z, i.e. the non-resonance condition l.(b), 
thus a similar expansion also yields Wk,i belongs in R). 

Now wf. 2 are solutions of the following system for fc = 1, • • • , N 
A?K,2)" + AIK.2) = m-s)+egl{s;w^^), 

iwym = o, K2)'(o) = o, 

with ^ wi + w^, w^^ {■■■ ,wl2, ■■■) and g%{s; ■■■ ,Wk,---)^ gl{s; ■ ■ ■ ,Wk,i + Wk, ■■■)■ 
The end of the proof of Lemma 16.41 is a straightforward generalization of the the proof of 

N 

Lemma 10 with the energy: 2E{wi,--- , wat) = ^ ((Ai)^(wfe)^ -|- (Afe)^Wfc) . □ 

For systems, we also have to work with linear combination of periodic functions with 
different periods and nonlinear function of such sum. So we work with the adherence in 
L°°{R,C) of spanje*'*'*, A € K}, namely the set of almost periodic functions C"p(M, C), and 
the Hilbert space of almost-periodic function is Llp{R,C), see with the scalar product 



1 '-^ 



{u,v) = linr - / u{t)v{t)dt. 

We give an useful Lemma about the spectrum of |u| for u G C"p(M,R). Let us recall 
definitions for the Fourier coefficients of u associated to frequency A: c\[u] and its spectrum: 
Sp [u] , 



[u] = {u,e'^') = lim u{t)e-'^*dt, Sp [u] = {\ e^, cx[u] ^ Q). (71) 



Lemma 6.5 [Property of the spectrum of \u\ ] 

If u £ Cqp(M,M), and if u has got a finite spectrum: 

Sp [u] C {±Ai, • • • ,±AAr}, with (Ai, • • • ,\n) being Z-independent, such that ^ Sp [u], 
then Xk ^ Sp [ \u\ ] for all k. 

Proof : The result is quite obvious for u^. We first prove the result for /(u^) where / is 
smooth. Then, we conclude by approximating \u\ by a smooth sequence /„(u^) = + , 

and using the L°° stability of the spectrum. 

Let E be the set of all Z linear combinations of elements of ^2 = {0, ±A^^ , k,j — 1, - ■ ■ N}, 
where 

^% = AfeibAj. 

Thus Sp [/(u^)] is a subset of E since Sp [u^] C S'2. 

Notice that X% = ±X% , A^, = 0, A+, = 2A, = A+ + A^^. . 
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Choosing /c = 1 for instance, so Ai ^ 0, it suffices to prove that Ai ^ E. 

Assume the converse, i.e., Ai G E. Then, for k < j, there exists some integers {c^j)k<j) such 

that: 

k<j 

Therefore, defining c^j, by ic^^ for k < j, we have: 



Using the Z-independence, with dkj = c^^- + c^T^ for k ^ j and rf^fc = 0, we have following 

system: 1 = Di — ^ dij, = 1?^ = ^ d^j , for all fc > 1. 

j j 

1 

Summing up, the iV — 1 last equations in — , and using the fact: djk = dkj modulo 2, we 

N N N 

have: = ^ Dk = ^ + 2 ^ d^j = ^ dij, then fi = 0, i.e. Di is even. It's impossible 

k=2 j=2 k<] 3=2 

since Di = 1. So Xi ^ E and the proof is complete. □ 



6.2 Bounds for the remainders 

Now, we prove each asymptotic expansion given in previous sections, i.e. we bound each 
remainders with energy estimates up to a large time. 

Proof of Proposition 12.11 : First we give the outline of the proof. 

Notice that all these computations only involve the function cos. Then, the only way to have 
a secular term in equations defining vi and V2 is a cos(s) in the right-hand side. So, the good 
choice of ai and a2, is enough to remove secular term with cos(s). Now, it suffices to control 

for large time. A computation shows that the remainder i?^ of equation (jl3p satisfies: 
\R,{s)\ < Ce{l + \r,{s)\) + \Xe\ivo,Vi + er,). 

Then, is like in Lemma 16. 3[ and the term comes from which is estimated by 
Lemmas 13.11 13.21 

More precisely, an exact computation of in equation (fTS]) leads to 



Re = Xsivo,vi + er^) + eH{vQ)re + ealv^ 



where af is a real constant, bounded uniformly for all e G [0, 1] such that 
(we)^ = ao + eai + £^ci2 + s^ct'^- From ((23l) we also have 

_ {vq + evi + £^rg)+ - [ivo)+ + H{vo){evi + e^r^)] 

£ 

_ ivo+evi) + -[{vo) + +eH{vo)vi] _^ {vo+evi+e^r^) + -{vo+evi)+ _ ^fj(^y^'jj.^ 

= Xs{vo, vi) - eH(vo)re + ege{s,r^), 

since u — > is 1-Lipschitz 

^, . {vq + evi + e'^re)+ - {vo + evi) + 

£ 

\9e{s,r,)\ < \r,\. 
So, with V;; — Vq + evi + E^r^, we can rewrite Rg as follow 

Re = Xe{vo,vi) + ege{s,re) + ealv'^. 
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Now, we can rewrite equation (jl3p in the following way 

-ao«+r,) = S{s) + f,{s)+eg,is,r,), 
S = a2VQ + aiv'l + H{vo)vi, 

fe = Xe{vQ,Vi) + £al{vQ' +£Vi'), 
9e = ge + S^al'Te, 

which allow us to conclude with Lemma 16.31 □ 



The proof for other propositions 12. 2[ 12.31 in section [2] are similar. 
We now complete the proof for the asymptotic expansions for systems given in section |4l 

Proof of Theorem 14.11 : As in the proof of Proposition 12.11 the same technique of proof 
is used component by component for Theorems 14.11 with similar energy estimates we can 
conclude with the Lemma 16.41 for system to control all r| . 
More precisely a complete computation of the remainder gives us: 

Ri = XeiakiVi - bk,Y^ akj[vij + er|]) + eH{akiv1 - bk) ^ afe^-r^ + ea|(t)|)", 

j j 

with notation of the proof of Proposition l2.1l and vf, — v''l + evl + £^r|. Let u, v, w be three 
functions, as previously, we have: 

Xe{u,v + £w) + eH{u)w — Xel^^i ''^) + + + ^ + 

and, since w — s- w+ is 1-Lipschitz: 

\e^^{{u + ev + e'^w)+ - {u + ev)+)\ < e\w\. 

Now, we can rewrite i?| as follow: 

j 

where g| is defined by 

j 

= Gkivl -bk+eY^ akjvlj , 
j 

Notice that g| satisfies |5|(s,rf, • • • ,r%)\ < ^ |afej||r||. 

j 

A key ingredient is the energy 2E = ^(ao(?'J,)^ + A^.r^) for the homogeneous system: 

k 

Lkfk — 0, k = 1, ■ ■ ■ , N and the for the inhomogeneous system: 

-LkVk = Sl{s) + f^{s)+sgkis,rl,--- ,r%), 

for = 1, • • • ,N, with 

N 

Sk = H{akiv1 - bkj^akjvj + a2V°" + aiVk' , 
j=i 

and ai, a2 are well chosen to avoid secular term when k = 1. Thus, all Sk are 2tt periodic. 
5*1 is not resonant with Li. The Xk are Z independent. We can apply Lemma 16.41 which is 
enough to conclude the proof. □ 
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Proof of Theorem 14. 31 : The proof follows two steps. First the solution for linear equations 
(j66p are bounded by o{t). Second, energy estimates are used to bound r^. 
At the end we prove remark [42l 

Notice that we do not use Lemmas 13.11 13.21 Indeed, we have no term with Xe- We only 
use that functions and uj? are Lipschitz, the Lemma 16.51 to identify resonant terms when 
6fe = and an energy estimate. But, since all modes are excited, the accuracy is weaker than 
the precision obtained in Theorem 14. li as in |24j . 



Step 1: the N problems ([66]) involves decoupled equations rewritten as follow with cj > 0, 

There is no resonance since ±w are not in the spectrum of S. But, Sp[S] is dense in R. 
Indeed Ai, - • • ,Ajv are Z independent. In general, we cannot expect that (p is bounded on 
the real line, see [B], but (p is less than 0{s) for large time. We can compute explicitly (f> 

0(s) — A cos {uj s) + B sin{ujs) + ijjls), 
Loip{s) = I S{(T)sa\{uo{s — <7))da 



sin(a;s) / S'(cr) cos(wCT)dcr — cos(ws) / S'(cr) s\n{LLi(T)da 



The condition ±aj ^ Sp[S] is lim s ^ I S{a) exp{±iu}a)da ~ 0. That is to say 

/ S{<t) exp{ztiuja)da = o(s) when s +oo, thus and cj) are negligible compared to s for 
Jo 

large time. 

Step 2: Let us decompose the remainder in the following way fl. = (pk + w^. From equation 



(j65p and the previous bound for (pk we have in variable t instead of Sk for convenience 

LkwUt) = 0{et) + {0{e^k) + 0{e\w'\))^0{et) + 0{e\w'\),) 

since (j>k{t) = o{t). Now, we remove the first part of the right hand side with w| = w| + Zj. 
and w| is solution of Lkwf. = 0{et). Classical energy estimates (or explicit computations as 
for (f>) yields to w'j,{t) — 0{et'^). Thus there exists a constant Ci > such that z| satisfies 

\L,zl\ < CieH' + e\z^\). 

Multiplying each inequality by summing up with respect to fc, integrating on [0,r], 

by Cauchy-Schwarz inequality, with D = 2C(min(Afc) + min(A;;)^) we get 

N 



E(T) = E(A?((z|)r + A?.(zD') 
fe=i 

< 2Ce2r2 r \^\[zl)'\{t)dt + 2Ce C \{z')' ■ z^dt 
Jo , Jo 



k=\ 



< De^T^'^^ E{t)dt + DeJ^ E{t)dt. 

Let Y{T) be / E{t)dt, thus Y{0) = and for aU t G [O, T], 
Jo 



E{t)^Y'{t) < De^T^-^y/Y{¥)+DeY{t). 
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Since [ = 2 In f 1 + wc obtain JWr) < eT^"" exp(DeT) and then 

Jo ^Vy + y \ A J 

E{T) < 2De^T'''exp{DeT). 

Finally r| = 0fc + w| + 2| = o(T) + 0{eT^) + 0{e^-^T^-^ exp{DeT)), so for any T, = o{e-^) 
we have in W^'°°{0, T,) for all T <T^ 

£r,{T) = o{eT,) + 0{£^T^) + Oie^-^T^-"^), 

which is enough to have the convergence in W^^^°°(0, Te). Furthermore r| satisfies the second 
order differential equation ((64|l which is enough to get the convergence in W'^'°° . 

About remarkl4.2l From Thcorcm l4.3[ this result its obvious. Let us explain why wc cannot 
go further up to the order e^. 

Unfortunately 5*^ is not periodic since is quasi-periodic for j ^ 1. Indeed, the following 
initial conditions 

4(0) = 0, K)'(o) = o, fc^i, 

yields to a quasi-periodic function, sum of two periodic functions with different periods 27r 
and 27rAi / Xk , thus a globally bounded function 



vlis) = ^i(s)- 01(0) cos (^^.) 



So we cannot apply Lemma [ 
Let us decompose Sk = Pk + Qk for k ^ 2 where Pk is periodic and Qk is almost-periodic 

Qk{s) = -H{akiVi{s) - bk)^akj<l>] {0} cos l^^s 



Let Wk be a solution of —LkWk = Qk then Sp[wk] G (J i"'"^ ^ ^ f' spectrum of Wk 

is discrete and there is resonance in the — 1 equations, —Lkrf. ^ Sk + ■ ■ ■ , k ^ 1 and the 
expansion does not still valid for time of the order e~^. □ 
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